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Abstract 

This article treats chaotic scattering with three degrees of freedom, 
where one of them is open and the other two are closed, as a first step to- 
ward a more general understanding of chaotic scattering in higher dimen- 
sions. Despite of the strong restrictions it breaks the essential simplicity 
implicit in any two-dimensional time-independent scattering problem. In- 
troducing the third degree of freedom by breaking a continuous symmetry, 
we first explore the topological structure of the homoclinic/heteroclinic 
tangle and the structures in the scattering functions. Then we work out 
implications of these structures for the doubly differential cross section. 
The most prominent structures in the cross section are rainbow singular- 
ities. They form a fractal pattern which reflects the fractal structure of 
the chaotic invariant set. This allows to determine structures in the cross 
section from the invariant set and conversely, to obtain information about 
the topology of the invariant set from the cross section. The latter is a 
contribution to the inverse scattering problem for chaotic systems. 

PACS Numbers 05 45 



1 Introduction 



Chaotic scattering with two degrees of freedom is fairly well understood [TJ 
[HIIU 02 [6] both for hard chaos (hyperbolicity) and soft chaos in terms of Smale's 
horseshoe construct [7J [8] . In a time- independent Hamiltonian system with two 
degrees of freedom, Smale's horseshoe displays the invariant manifolds which 
qualitatively separate the dynamics. On an appropriate surface of section these 
invariant sets become smooth curves. This has two advantages:!) we can study 
objects easily if they can be drawn on paper, and 2) they separate the phase 
space dynamics if they close or go to infinity in both directions. 

Scattering functions jH O HO] and cross sections [TTJ [T^] have been analysed 
both statistically and geometrically, and the inverse scattering problem has 
been tackled [HJ H3] (for general background information see the text books 



Now the effort must be directed towards higher dimensional systems, as they 
are relevant in astrophysics (TTl [18] or chemistry (T9J [20J, [2T] . How do we proceed 
to generalise this success in the description of chaotic scattering to more than 
two degrees of freedom? There are previous important steps particularly by 
Wiggins [2"2 l l2"5 | l2l], but also by Ott 25J. The former develops a formal theory 
that is difficult to apply but gives important foundations of the problem. The 
latter show, that a straight forward generalisation of the three-disc system in 
triangular configuration to a four sphere problem in tetrahedral configuration is 
feasible, but the invariant manifolds are low dimensional and thus not relevant 
so some physical problems, as they are difficult to detect. We have suggested in 
an earlier paper that it can be useful to have a gradual approach to a change in 
dimension by breaking a continuous symmetry of the problem [26 . 




In the present paper we shall follow this last suggestion and apply it to the 
simplest possible three degrees of freedom configuration, with two bound and 
only one open degree of freedom. Despite of this strong restriction, we are now 
definitely beyond the situation, where Smale's theory can be used, and systems 
of that kind are not without interest. We might consider guiding channels with 
one open degree of freedom, or we can consider equivalently three interacting 
particles in one dimension, under conditions where asymptotically one pair must 
always be bound. 

We shall focus on the topological structure of the invariant set rather than 
on details of trajectories, because the former properties tend to be generic , 
i.e. robust and embody the most relevant information about the chaotic system 
[T5j . We are thus looking for some generalisation of the horseshoe to higher 
dimensions. A simple minded generalisation does not exist, as can be seen 
from the work of Wiggins [53] . The basic reason is that the dimension of the 
invariants manifolds of hyperbolic points is too low. 

We proceed by analysing the transition from a system, which has an 0(2) 
symmetry group in the three dimensional position space, reducing the problem 
for a fixed value of the angular momentum to a system with two degrees of 
freedom. We shall then introduce in the reaction region (i.e. in the region 
where the asymptotic Hamiltonian is not valid) a symmetry breaking term. 
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This will allow us to start our analysis with a continuous stack of horseshoes 
corresponding to different values of the angular momentum and observe how the 
structure evolves as the invariance breaks down. We then study the implications 
for scattering functions and for the doubly differential cross sections. The fractal 
structure of singularities in the latter, i. e. the rainbow structures, is directly 
related to the structure of the chaotic set, leading to an additional important 
contribution to the corresponding inverse scattering problem. Indeed we shall 
also see, that the basic idea developed is not limited to three degrees of freedom 
nor to one open degree of freedom, but these aspects will only be commented. 

We shall illustrate our procedures with three examples. First a channel 
with harmonic walls in two directions that deforms to some more complicated 
potential in some compact region of configuration space. This is equivalent to 
two pairs of particles bound by harmonic forces and interacting between each 
other weakly. The second is a bottle shaped billiard that is either connected 
to a channel billiard or has an opening that separates the interior from the 
exterior in some plane in configuration space. Finally for numerical convenience 
we shall use a 5 kicked system of two degrees of freedom, which is topologically 
equivalent to the Poincare maps of the first system. 

As far as the physical implications are concerned, it does not seem easy 
to emulate systems with such restrictions and no significant friction, though 
it might become feasible with Bose Einstein condensates at some point in the 
future. The alternatives in celestial mechanics imply higher dimensions and/or 
more complicated asymptotes, as is the case with the breaking of the Jacobi 
invariant in ref [26] ■ Yet the results are certainly useful for semi-classical con- 
siderations, such as e.g. the ones used in |27j to design and interpret a microwave 
experiment. 

The paper is structured as follows: in the next section we shall define the 
class of three degree of freedom systems, which we shall consider as well as the 
specific examples we shall use. Next we proceed to construct the invariant sets 
and study them as a function of the symmetry breaking. In section four we 
derive the scattering functions, that is the behaviour of outgoing asymptotes as 
a function of the incoming ones. This will allow us to discuss the rainbow effects 
in the doubly differential cross sections in section five. We finally proceed to 
comment on the scope of our results as well as implications for semi-classics. 

2 The class of systems considered and the model 
systems 

We consider a scattering system with three degrees of freedom where one of 
them is open and the other two are closed. We imagine that the first degree of 
freedom is a translational degree of freedom between projectile and target, the 
second degree of freedom is a vibrational degree of freedom and the third one is 
a rotational degree of freedom. 

The first particular model system of this class which we use is a rotationally 
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symmetric channel containing an additional short range potential representing 
some obstacle (the target) sitting in the otherwise empty channel. Accordingly 
the total Hamiltonian of the system splits, as always in scattering systems, in a 
free Hamiltonian Hq describing the asymptotic free motion ( in this particular 
case it is the motion in the empty channel), and an additional interaction W 
which is the short range scattering potential, i.e. 



H = H + W (1) 

We imagine that the channel runs straight in one direction and we use the 
coordinate q along the channel, p is the momentum conjugate to q. In the 
transverse direction the motion of the projectile is confined by the channel 
and we assume that the potential which represents the walls of the channel is 
quadratic in the transverse coordinates. Because of the rotational symmetry 
of the empty channel one natural choice is to use cylindrical coordinates p and 
9 and their conjugate momenta p p and L and also to include into the free 
Hamiltonian the quadratic potential and the transverse kinetic energy in these 
coordinates. Since the radial motion is an oscillation, it is sometimes more 
convenient to use action and angle variables I and </> for the radial degree of 
freedom. We will switch freely between these two possibilities according to which 
one is more convenient at the moment. The asymptotic Hamiltonian written in 
cylindrical position and momentum coordinates results as 

H =p 2 /2+p 2 p /2 + L 2 /(2p 2 ) + p 2 /2 (2) 

We choose units of time such that the transverse oscillation frequency in the 
empty channel is one. To describe the obstacle we use later as model of demon- 
stration the potential 

W^-^tB (3) 

where 

D 2 = q 2 + p 2 (sm 2 9 + (l + A) 2 cos 2 6) + 1 (4) 

The extra constant, taken as one, avoids singularities in the potential. Very 
important for the following is the parameter A which measures the deviation 
of the interaction potential from rotational symmetry. For A — we have 
perfect symmetry and since the free Hamiltonian Ho is also symmetric, the total 
system is symmetric, the angular momentum L is conserved and the system can 
be reduced to two degrees of freedom. For the value L = the symmetric 
system reduces to the one used in This reference also contains a physical 
interpretation of similar potential models. 

Next we need labels for the asymptotes of the system, i.e. the trajectories 
of the free motion described by Hq. The optimal possibility is to use 5 indepen- 
dent quantities which all are constant along the asymptotic trajectories. Since in 
any Hamiltonian system the energy E is conserved, both under the free asymp- 
totic dynamics and under the full dynamics, we use E as one of these labels. 
Asymptotically the translational degree of freedom decouples from the other 
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ones and the momentum p becomes constant. We use p as second asymptotic 
variable. The sign of p indicates into which direction the asymptotic motion 
runs and therefore a comparison between the signs of initial and final p distin- 
guishes transmission and reflection. In addition, since H is independent of 9, 
the conjugate variable L is conserved under the asymptotic dynamics and it is 
convenient to use it as third label for asymptotic trajectories. To distinguish 
the various asymptotic trajectories with the same values of E, p and L we need 
two additional labels giving the relative phase shifts between the translational 
motion along the channel and the other two motions. The systematic choice 
for these relative angles are the reduced phases constructed along the following 
idea. In action variables we find for our particular system 

H (p,I,L)=p 2 /2 + 2I+\L\ (5) 

Then the asymptotic equations of motion for ip and 9 are 

d(j) dH a 

Tt = ^r =w ^ L) = 2 (6) 



d9 dH 

- = _ = a*(I,L)=±l ( 7) 

The sign appearing in equation [7] is the sign of the angular momentum. We 
define the reduced angles ip and \ belonging to 9 and tf> respectively as 

^ = 6-w e (I,L)q/p = 6Tq/p (8) 

X = 4>-u 4 ,(I,L)q/ P = 4>-2q/p (9) 

The sign in equation [8] is minus the sign of the angular momentum. A short cal- 
culation shows that these two reduced angles are constant under the asymptotic 
dynamics. 

An equivalent possibility is the following: We use the values of the two 
coordinates <f> and 9 at the moment when the absolute value of q reaches some 
very large value (in the numerical examples, we use ||<7|| = 8). We call these 
two particular values which serve as asymptotic labels x an( i "0 again. We use 
the index in for initial asymptotic conditions and the index out for the labels 
of outgoing asymptotes. 

Since it is simpler to handle iterated maps instead of flows, we will use 
Poincare maps to represent the dynamics of the system and to explain many 
ideas (an elementary explanation of the concept of Poincare maps is found 
in section 2.5 of |28j). For the channel problem an appropriate intersection 
condition for the Poincare map is to take the maximum of the cylindrical radial 
coordinate p. This choice has the advantage to coincide with the choice made 
for a billiard model which we also use later. Almost all trajectories intersect this 
surface transversely an infinite number of times and asymptotically the return 
time to this surface becomes constant, note the constant value of w^(J, L) = 2 
in equation [6] The only exceptions are trajectories with action 1 = which 
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run along spirals of constant value of p without any radial oscillations. In the 
domain of the map we use the canonical coordinates q, p, 9, L. 

An alternative version of the reduced angles ip and X f° r the map is obtained 
as follows: For trajectories with the value p of the asymptotic longitudinal 
momentum select an interval of large absolute values of say the interval 
[Q, Q + \p\] where Q is sufficiently large to guarantee that the trajectory is 
already in the asymptotic region. Observe the trajectory of the map until it 
steps into this interval. For this point define x = 27r(|g| — Q)/\p\ and for i/j 
take the actual value of in this point. Note that the initial point and the final 
point of the selected q interval can be identified for the purpose of asymptotic 
labelling since they describe the same trajectory (one is the image of the other 
point under the action of the map) and lead to the same value of the angle 
X since multiples of 2tt are irrelevant. Therefore this line of initial conditions 
has the topology of a circle. An analogous construction also works well for the 
Poincare maps of other model systems. 

The second model of demonstration used is a billiard model, it describes 
scattering in a bottle. The boundary of the bottle will be defined with the aid 
of the following functions: 



f( s ,1-Wl-q*, forg<0 

fW = 1 '5 , 2,1, + u • ( 10 ) 

aoq + aiq + 1.4, otherwise. 



r(q,6) = 1 + Acos(0)cos(^) 2 (11) 
The boundary itself is given by 

p( q ,e) = r(q,9)f( q ) (12) 

Two dimensional billiards with the border described by function [10] have 
been studied in [22] and a similar model in [3D] . 

The values of the constants a; are chosen to give a smooth boundary and 
to lead to an unstable quasi periodic orbit in the plane q = so that there 
is a complete binary horseshoe for the reduced rotationally symmetric case for 
L = 0. Convenient values turn out to be: 



ao = ^(3.49)" 1 / 4 (13) 
«i = -0.7 (14) 
go = (3.49) 1 / 2 (15) 

Also this model has a parameter A which gives the distance from rotational 
symmetry and again for A = we have perfect symmetry and the system can 
be reduced to one with two degrees of freedom. 

The billiard dynamics is the usual one with elastic reflection on the wall and 
free motion in the inside. The energy is a constant and proportional to the 
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square of the momentum. Therefore we can set the velocity ||u|| = 1 without 
loss of generality. The Poincare map of the bottle model is the usual Birkhoff 
map on the wall of the billiard (these are explained in [31]). We only use weak 
deformations of the bottle such that this map coincides with the intersection 
condition of maximal cylinder radius p. Accordingly we use in the domain of 
the map the same coordinates q, p, 9 and L as in the channel model. 

The Birkhoff-Poincare map of this system has a binary horseshoe instead of 
the ternary one of the other two examples, and it has a surface of no return at 
a finite value of q, which makes is qualitatively different. Nevertheless we will 
see that it belongs to the same class of systems. 

Since it is a lot more convenient and faster to investigate maps instead of 
flows we construct as third model a closed form example of a map acting on the 
4 dimensional domain with coordinates q, p, 9 and L. This model is constructed 
in a way that it can serve as a prototype Poincare map for the whole class of 
systems considered in this article. It is based on the usual scheme of kick and 
free flight. The generating function for the kick is: 

G(q, 9,p, L) = qp + 9L+ (L max - L)(l + Acos9)V(q), (16) 
with potential function 

V(q) = -e-* 2 . (17) 
We define the force function on the q coordinate as 

F(q) = -dV/dq = -2qe- q \ (18) 

And the kick map is implicitly given by: 



- dG dG 

?= 7* P =Dq (19) 

9 ~dL L -m- (20) 

To construct a complete step of the map, the particles will perform first half 
of the free flight: 



q' = q„+Pn/2 (21) 
9' = 6 n + L n /2 (22) 

p' = Pn (23) 

L' = L n (24) 

Afterwards, we give this auxiliary coordinates the kick, given explicitly by 
the next set of equations: 
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q" = q' (25) 

6" =9' -(l + Acos6')V(q') (26) 

P =P + (L max - L ) 1 + AV{ql)s - m6 , (27) 

AF(?)sin(9' + 1 1 j 
And finally half a free flight is applied again: 



q n+1 = q" + p"/2 (29) 

0„+i = 9" + L"/2 (30) 

Pn+i=p" (31) 

= L" (32) 

This concludes the action of the map. Every step is symplectic, and the map 
behaves qualitatively as the Poincare map of the first example. 

We use mainly the map model for the presentation and explanation of our 
ideas. Below we compare some results of this map with the results of the other 
two models in order to convince ourselves that the map results are representa- 
tive. 

3 The topological structure of the chaotic in- 
variant set 



Since the interaction potential of equation [3] or of equation 17 is negative, there 



are trajectories with negative energy yet arbitrarily close to zero, going out ex- 
tremely far and returning. Therefore our first example system has no points 
without return. Its outermost periodic orbits lie at q — ±oo. They are trajec- 
tories oscillating transverse to the channel and rotating at constant value of 
where q is arbitrarily far away. Since any displacement of q leads to an equivalent 
transverse oscillation, these transverse trajectories come in a whole continuum 
of copies. In the Poincare map they lead to fixed points at q — ±oo, p = 
and these points are parabolic, i.e. neutrally stable in linear approximation. 
However the non-linearities of the map at finite values of q make these points 
non linearly unstable and they have stable and unstable invariant manifolds. 
The stable manifolds consist of trajectories which go out to infinity monoton- 
ically, i.e. the absolute value of q increases monotonically while at the same 
time the value of the momentum p converges to zero. Asymptotically all energy 
of the motion goes to the transverse motion. The unstable manifolds consist 
of trajectories doing the same under the time reversed motion. In this sense 
the trajectories belonging to the invariant manifolds of the points at infinity 
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converge to the periodic trajectories described above. In the domain of the map 
these manifolds are the separatrix lines between motion going monotonically 
away and motion which returns. 

The map model described as third example of demonstration in the previous 
section is in this respect similar to the channel. It also does not have values of 
q of no return. 

In contrast the bottle model has a surface of no return, the surface of the 
bottle neck. The trajectories lying on this surface forever serve as outermost 
fixed points in the Poincare map. 

The set of trajectories staying in the plane of the bottle neck and also the 
trajectories with p = in the asymptotic region of the channel form an example 
of what Wiggins 23] calls a normally hyperbolic invariant manifold, abbreviated 
NHIM. Let us explain this set first for the case of the bottle, since in this case 
we have a NHIM in its original form. 

For fixed energy, all the orbits which stay forever in the plane of the bottle 
neck form a two dimensional continuum. First, there are such orbits for all 
possible values of the angular momentum and second to each such orbits there 
also exist corresponding orbits rotated by an arbitrary angle. These trajectories 
are neutrally stable under perturbations of initial conditions which keep them in 
the bottle neck plane. In the domain of the Poincare map all these trajectories 
form a 2 dimensional surface, which is the NHIM in the domain of the map. In 
the direction perpendicular to the bottle neck all these trajectories are unsta- 
ble and they have stable and unstable manifolds. The union of these invariant 
manifolds taken over the whole 2 dimensional continuum of bottle neck trajec- 
tories form the stable and unstable manifolds of the NHIM which we will call 
W s and W u . They are 3 dimensional surfaces in the 4 dimensional domain of 
the map. Therefore these surfaces are dividing surfaces, they divide trajectories 
which pass the bottle neck from trajectories which return, for more details and 
consequences for scattering see [22], or for applications see [24] . 

In the example of the channel we have a slightly modified version of a NHIM. 
Trajectories in the empty channel or in the asymptotic region of the channel 
having p — stay forever in the plane of a fixed value of q, they are trajectories 
in the two dimensional oscillator potential forming the empty channel. In a 
harmonic oscillator all trajectories are periodic, they are ellipses. The ones for 
a fixed value of the total energy are a two dimensional continuum where each 
individual trajectory can be distinguished by its value of the angular momentum 
and by the orientation of the ellipse. Again such trajectories are neutrally stable 
under perturbations which keep the value p = 0. The difference between the 
channel and the bottle example is that in the channel case the trajectories 
which form the NHIM are also neutrally stable in q direction, at least in linear 
approximation. If we include the non linearities caused by the asymptotic tail 
of the attractive scattering potential, they become unstable and have stable 
and unstable manifolds. The union of these invariant manifolds forms again 
the invariant manifolds of the whole NHIM. They are dividing surfaces which 
separate trajectories going out monotonically from returning trajectories. In this 
respect the map model behaves like the Poincare map of the channel model. 
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To investigate the chaotic set of the models we start with the case of A = 
and later on break the symmetry. For A = the 8 dependence of the dynamics 
decouples from the rest of the dynamics, L becomes a conserved quantity and 
the map reduces to a 2 dimensional map in the coordinates q and p, where L 
serves as a parameter. Because of its importance in what follows we show this 
reduced map in explicit form for the map model with the potential function 



from equation 17 



q n +i = Qn +Pn - (q n +p„/2)exp(-(g„ + p n /2) 2 )(L max - L) (33) 
p n+1 =p n - 2{q n +p„/2)exp(-(g n + p„/2) 2 )(L max - L) (34) 

(35) 

The L value determines the strength of the force of the kick, becoming 
negligible for L sa L max . This corresponds to all energy spent on rotational 
motion, and almost no development for the horseshoe, as can be seen on the 
sequence of plots of the homoclinic tangle for this reduced map in figs, [I] 

To construct the horseshoe we plot the stable manifold and the unstable 
manifold of both the fixed point at +oo and of the one at — oo. To plot the 
whole interval q £ [—oo, +oo] on a finite range we use the horizontal coordinate 
z = tanh(g) which is convenient to get informative plots (figure [2]). The inter- 
section points between a stable and an unstable manifold are trajectories which 
converge forward in time and also backward in time to a fixed point of the map. 
An intersection between the manifolds of the same fixed point is called homo- 
clinic point, and an intersection between the manifolds of different fixed points 
is called heteroclinic point. The whole structure of the manifolds is called ho- 
moclinic/heteroclinic tangle, many times only the expression homoclinic tangle 
or simply tangle is used for short even when it also contains heteroclinic points. 
The existence of a homoclinic tangle is the topological criterion for chaos in 
the system. It implies the existence of an uncountable set of unstable trajecto- 
ries. For more details on the role of homoclinic tangles and the importance for 
scattering and transport problems see the book by Wiggins [2"3"] . 

Because of symmetry in our example the manifolds of point — oo are obtained 
from the ones of the point at +oo by inversion about the origin. The time 
reversal symmetry in the maps permits also to obtain the stable manifolds from 
the unstable ones by the reflection p — > —p. From those symmetry properties 
also follows the existence of a trajectory oscillating transverse to the channel at 
q = 0. It leads to a fixed point of the map at q = 0, p = 0. We call this fixed 
point the inner fixed point. In total the map has three fixed points where the 
outer two are symmetry related and therefore we call the resulting homoclinic 
structure a ternary symmetric horseshoe. 

When changing L we have the usual development scenario starting from a 
complete horseshoe for L — up to a parabolic line when L reaches its maximal 
limiting value L max allowed by total energy. In our particular example we have 

(36) 
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(a) L = 0.00 





Figure 1: The Poincare Maps for the discrete dynamical system in the axial 
symmetric case (A = 0). The L parameter regulates the degree of development 
of the horseshoe. Notice also reduction of the phase space occupied by the 
horseshoe as L increases. The line C on 1(a) it's choosen so that it's preimage 
is an adecuate domain for scattering functions. 




-0.5 

(a) L = 0.00 




-0.5 
(b) L = 2.60 




-0.5 
(c) L = 5.71 



Figure 2: Same as fig. [T] but after compactifying the q axis using the function 
arctan(</). 12 



Next we define the fundamental area R for the horseshoe. Let us start with 
the local branches of the manifolds of the two outer fixed points and let them 
grow longer continuously in a completely symmetric pattern. Then at some 
length we obtain the first intersection points, for symmetry reasons they are 
located at q = 0. At this moment let us stop the process of growth. The result 
is a curvilinear quadrilateral which we call the fundamental rectangle R. 

The intersection pattern between some local segment of an unstable man- 
ifold (e.g. the piece between an outer fixed point and the next corner of R) 
with the whole fractal bundle of stable manifolds already contains the com- 
plete information of the whole tangle formed by the manifolds of the outer fixed 
points, and the knowledge of such local intersection patterns is sufficient. We 
will make use of this idea for the reconstruction of important properties of the 
tangle from scattering data. The use of these for determining the development 
of the horseshoe is profusely explained in [13] . 

In figur^l] a) we included an additional thick line labelled C which runs 
parallel to one of the local segments of unstable manifolds but just outside of 
R. Note that the intersection pattern between C and the bundle of segments 
of the stable manifolds coincides exactly with the intersection pattern between 
the local segment of the unstable manifold and the bundle of segments of stable 
manifolds. Accordingly also the intersection pattern along line C contains all 
important information about the whole tangle. The iterated preimages of the 
line C play an important role in the following as asymptotic initial conditions 
for scattering trajectories. 

Now we will discuss the chaotic structure of the three degrees of freedom 
system. We begin with A = 0. Then the surfaces of the various values of L arc 
dynamically independent and we get a fractal structure in the full dimensional 
domain of the map with its coordinates q, p, L and 9 in a two step process. First 
we pile up the tangles for the various values of L obtaining a fractal tangle in a 
three dimensional embedding space and second we form the Cartesian product 
of this object with a circle representing the fourth coordinate 9. 

For the general three degrees of freedom case we need an argument of ro- 
bustness of this structure under moderate deformations and moderate breaking 
of the rotational symmetry. This argument follows from generic transversality 
properties. 

In a system of the form described by the equations [T] and [2j or in the map 
we obtain a family of horseshoes as depicted in figure [T] in the symmetric case. 
The development of the horseshoe is indicated by a parameter 7 which measures 
the degree to which the horseshoe is developed. 7 = is the parabolic line and 
7 = 1 is the complete ternary symmetric horseshoe, for more explanations and 
schematic pictures see |13j . In our present case 7 is a monotonic function of L. 
Accordingly also the value of L orders the pile of two dimensional horseshoes 
naturally from the complete case, with L near zero, to the integrable system, 
with L near L max . When we view this family of systems as just one geometrical 
object then the family of invariant manifolds form smooth two dimensional 
surfaces in the three dimensional embedding space with coordinates q, p and L. 
For short we call this structure the three dimensional tangle. Now let us restrict 
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Figure 3: A plot of the intersection of stable and unstable manifolds in the 
(p, q, L) space when A = 0. The L value gives the degree of development of the 
horseshoe for every (q,p) plane. After stacking these in order, we obtain two 
intersecting two dimensional manifolds. In the end we consider a direct product 
with a circle for the 9 coordinate. 



our attention to intersections between the unstable and stable manifolds. The 
basic pattern is a transverse intersection between two two dimensional surfaces 
in a three dimensional embedding space as shown schematically in figure [3j 
A point of interest is that the tangential intersections of the lower dimensional 
horseshoe become extremal points in L direction of the intersection curve shown 
in figure [3) The whole curve belongs to the "hyperbolic component" of the 
invariant set. 

So far we have ignored the 9 coordinate. To include it for the symmetric 
case of A = we simply form the Cartesian product of the structures described 
above with a circle representing this cyclic coordinate. Thereby the manifolds 
become three dimensional surfaces in the four dimensional embedding space. 
The elementary transverse intersection structures between the surfaces W u and 
W s are smooth two dimensional surfaces. Each one of their points has one 
dimensional unstable and stable directions and two neutral directions. Because 
the manifolds themselves are folded to fractal structures we obtain a fractal 
repetition of the elementary intersection structures which we call the tangle in 
the four dimensional domain of the Poincare map. The dimension of this tangle 
can be anything between two and four depending on how exactly the manifolds 
turn and fold. For an argument we present later, remember that this dimension 
is larger than two. 

Finally we consider the changes implied by a moderate breaking of the ro- 
tational symmetry. To understand this let us return to figure [3] and imagine 
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its product with a circle. Then we see the transverse intersection of two three 
dimensional surfaces in a four dimensional embedding space. Because the in- 
tersection is transverse it is structurally stable. This means, under small defor- 
mation the intersection pattern remains qualitatively the same. This argument 
applies to all the transverse homoclinic intersections appearing in the develop- 
ment scenario of the horseshoe. However, it does not apply to the non hyperbolic 
structures near the surface of KAM islands. 

So, when the parameter A is slightly different from zero, we can understand 
that relevant qualitative properties of this structure remain unchanged. The 
structure of the whole four dimensional tangle is therefore robust. This stability, 
together with the transversality of the intersection, is inherited from the leaves 
of the symmetric system in a way which is only possible for systems which can 
be connected continuously to the symmetric systems and which are not far from 
the symmetric case. This may appear restrictive at first, but even so we can 
cover a broad range of physically relevant problems. 



4 Scattering functions 

Scattering functions are interesting objects in their own right and we need them 
as auxiliary concepts to explain and motivate the ideas on the singularities in 
the cross sections, which are more accessible to experiments. 

A scattering function gives one of the outgoing asymptotic labels, or several 
of them, as function of initial asymptotic conditions. It is advantageous to 
study outgoing action like variables as function of initial angular variables, this 
is exactly the version of the scattering functions we need to understand the 
properties of the cross sections. As explained before, every asymptotic trajectory 
in the Poincare map can be labelled by p, L, tp and \. It is always understood 
that the total energy E is kept fixed at one particular value. 

The scattering function which we will now describe in detail gives the out- 
going momentum p out and the angular momentum transfer AL — L out — Li n as 
function of Xin an d tpin f° r fixed values of Pi n , L in and E. The domain of this 
function is the 2 dimensional torus with coordinates Xin and ipi n . Periodicity 
of the function follows thereby. 

First let us study the symmetric case A — and show some numerical 
examples for the model map. In this case the angle ip is irrelevant and because 
of angular momentum conservation the function AL is identically zero. Then 
the domain of the scattering function is an interval of length 2tt in the angle Xin 



and this interval is represented by the line C in figure 1(a) Remember that C 
is the iterated image of an asymptotic line of fixed p where q changes over an 
interval of length p which corresponds to an interval of Xin of length 2ir. 

As a consequence the scattering function has singularities on a fractal set 
and intervals of continuity in between. The singularities correspond the inter- 
section of the line of initial conditions with stable manifolds. Since preimages 
of intersections with invariant manifolds are again intersections with the same 
manifolds, the intersection pattern of this line of initial conditions coincides with 
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the fractal pattern of intersections of line C with the stable manifolds of the 
horseshoe. If the asymptotic part of the trajectory starts on a stable manifold 
of the horseshoe and the actual scattering trajectory converges to a periodic or- 
bit and does not come out of the interaction region with a longitudinal kinetic 
energy larger than zero, then the outgoing asymptotic conditions are ill defined, 
and therefore the scattering function has a singularity. The pattern of singular- 
ities is the pattern of intersections of the line C with stable manifolds and is a 
smooth image of the fractal pattern in the tangle. In this way asymptotically 
obtained scattering functions carry the information on the topological structure 
of the tangle sitting in the interaction region (T3] • 

We have to ensure that the line of initial conditions maps to a line equivalent 
to line C in the horseshoe plots and not to a line intersecting the outer tendril 
further out, where it does not scan the whole structure. As we move line C 
further away from R we loose more and more intersection points in the line, until 
it is so far out that it does no longer intersect the bundle of stable manifolds at 
all. Here it is essential to use a value of pi n sufficiently small. The actual line 
used for the numerical examples of the scattering functions and cross section 
data will be the 312th preimage of the line C shown in the figure, which is 
sufficiently far out to be considered asymptotic. This line is given by 

T- 312 C = {(q,p)\\q £ [-6.962, -6.912], Pm = 0.05}, (37) 

where T is the Poincare Map. If the line of initial conditions intersects the 
fractal structure only partially, then it still would contain the complete infor- 
mation because of the self-similarity of the fractal structure. However, then the 
reconstruction would pose some additional technical problems which we want 
to avoid. Typical structures for the scattering functions on a one dimensional 
domain are presented in figures [4] and [5] 

In a given plot of a scattering function we recognise intersections with the 
stable manifolds as initial conditions which lead to p ou t — 0. They represent 
boundaries between transmission (p ou t > 0) and reflection (p ou t < 0). 

In preparation for the asymmetric case we have to include the role of the 
coordinate ip. For A = the fractal of singularities in the two dimensional 
domain of the scattering function is the Cartesian product of a one dimensional 
fractal along the x direction described above with a circle in i\> direction. The 
intervals of continuity have the structure of strips which run in the ip direction 
around the domain. The natural domain of the higher dimensional scattering 
functions is the product of the preimage of the line C described above with the 
circle of tp values. It is a two dimensional torus if we recall that the initial and 
final point of the line C should be identified. 

Let us break the rotational symmetry and remember the robustness men- 
tioned at the end of the previous section. For A small the important intervals 
of continuity have still the same qualitative structure, they are only deformed 
continuously. Only the non hyperbolic dust around KAM islands is changed 
immediately for 4^0. In this sense the qualitative structure of the fractal of 
singularities is still close to a product of a one dimensional fractal with a circle. 
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Lj„ = 0.000 




Figure 4: The p out Scattering function for A = (axial symmetric case). The 
independent variable is qi n , the entrant? phase. 



k„ = 0.000 




-6.945 -6.94 -6.935 - 6.93 - 6.925 - 6.92 -6.915 



(a) 



Ii„ = 0.000 




-6.9409 -6.9409 - 6.94089 - 6.940S9 - 6.940SS -6.940SS -6.940S7 



(c) 

Figure 5: An illustration of the self similarity of the scattering function p out 
with parameter L in — 0.00. This case lj^s a complete horseshoe, so the function 
reveals the structure of a Cantor Set in the singular values. 



With increasing value of A intervals of increasing size and lower level in 
the fractal hierarchy are changed qualitatively in so far as they are disrupted 
into fragments which no longer have the structure of stripes winding around the 
domain in ip direction. Only for A « 0.5 the last large intervals of continuity are 
destroyed. Thereby the last remnants of the product structure of the fractal are 
lost and it is transformed into a truly higher dimensional fractal. In the present 
paper we shall not investigate this new structure, but restrict our consideration 
to smaller values of A. 

So far we have seen that all topological information on the chaotic tangle is 
contained in asymptotic scattering data, and that a measurement of scattering 
functions and their analysis is one way for the asymptotic observer to obtain 
this information. The scattering functions present a kind of shadow image of 
the chaotic invariant set cast into the outgoing asymptotic region. However 
the measurement of scattering functions is difficult or impossible in many real 
scattering experiments and therefore we now have to go one step further to cross 
sections which are the quantity measured in most of the traditional scattering 
experiments. 

5 The cross section 

The measurement of scattering functions needs the control over the canonically 
conjugate variables p, L and \i ^P- I n classical mechanics this could be done 
in principle. However, it is not done in the usual scattering experiment. More 
important, in quantum mechanics this preparation is forbidden even in princi- 
ple, and if we want to develop concepts and ideas which have some chance to be 
transferred to quantum dynamics, we can not use quantities needing the simul- 
taneous preparation of canonically conjugate variables. What is done in most 
scattering experiment is the following: One half of the phase space variables 
(normally the actions) is prepared as sharp as possible and the other conjugate 
half (the conjugate angles) is completely unspecified. In our case this means: 
For fixed total energy E we prepare p in and L in to specific values and do not 
have any control over the reduced angles Xin an d ipi n - Their values have a 
distribution with constant density over the whole domain, T 2 . The detector 
measures the final values p ou t and L out for each outgoing particle and we mon- 
itor the relative probability a(p out , AL) to find a given combination of values 
of Pout and AL, normalised by the incoming flux. This relative probability is 
called the doubly differential cross section. For more general information on 
cross sections see some text book on scattering theory as e.g. [33] or [33], for 
an application consult |12j . 

In the previous section we have seen that the scattering function with values 
Pout{Xin,ipin), AL{Xin, V'm) contains the fractal structure of the horseshoe. If 
we can measure this function directly then we have the necessary data to re- 
construct the pattern of the horseshoe and our version of the inverse problem 
is solved (13j . If we can only measure cross sections we must find out how the 
fractal pattern of the scattering function is transferred to some recognisable 
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pattern in the cross section. 

In the beam of incoming projectiles the weight on the Xim^in torus is con- 
stant and the scattering function maps this initial weight on some interval in 
the range of this function having the coordinates p out and AL and the image 
weight is exactly the cross section a(p out , AL). Therefore to get the value of 
a for a particular combination of values of p ou t and AL we first search for all 
preimages Xki ^k of the image point. Each preimage gives the contribution 

g k { PouU AL) = l/\Detd(p,AL)/d( X ,i>)\. (38) 

Then the value of the cross section is just the sum of these weights over all the 
contributing preimage points, i.e: 

a(p out ,AL) = ^ 9k(Pout, AL) (39) 



From equation 38 we see that the cross section has singularities at image 
points where the Jacobian of this map is zero, i.e. points which are locally not 
invertible. They are caustics of the projection of the graph of the function into 
the image space. The corresponding singularities in the cross section are called 
rainbow singularities, for their effect on light scattering off water drops. Now we 
shall focus on these singularities. Each interval of continuity of the scattering 
functions leads to one copy of a typical rainbow singularity in the cross section. 
This gives a possibility to see the fractal structure of the chaotic invariant set 
in the cross section. For systems with two degrees of freedom this idea has 
been worked out in detail in [TTJ [T^] and in this section we explain the higher 
dimensional generalisation. 

The determinant of the derivatives of the scattering function is zero exactly 
for such values of p ou t and AL for which two trajectories fall together and 
disappear. Accordingly these singularities are the lines across which the number 
of contributing trajectories, i.e. the number of preimages changes by 2. Now 
we will construct a simple analytical normal form for these singularity lines 
and compare it with numerical results for the three model systems presented in 
section 2. 

First we need an analytical model for the scattering function in one interval 
of continuity which is as simple as possible but still produces the generic form 
of singularity in the cross section. As explained in section 4, the domain of the 
scattering function is the 2 dimensional torus with coordinates Xin and ipin ■ For 
small deviation from symmetry a typical interval of continuity is a strip running 
around the torus in ip direction. In x direction the strips run in the symmetric 
case over a limited range only, let us say from xo ~ $ U P to xo + <^ where xo 
is the position of the middle of the strip. The value of p ou t is maximal in the 
middle of the interval, let us say it has the value po m the middle, and it goes 
to zero on the boundaries of the interval of continuity. Accordingly the simplest 
model function for p out in the symmetric case is 

Pout =Po~ (Xo -Xin? (40) 
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where the variables are scaled such that 8 2 = po- For the asymmetric case we 
add a tp dependent term. This term must be 2tt periodic and for simplicity we 
take the lowest Fourier contribution only, resulting in 

Pout =Po~ (Xo ~ Xin) 2 + bcos(ip m ) (41) 

for the general case with broken symmetry. Here for each given value of "0 the 
variable x 1S restricted to such values which result in a positive value for p ou t- 
This restriction defines the deformed interval of continuity for the case of broken 
symmetry. 

The function AL is identical zero in the symmetric case. For the case of 
broken symmetry we will consider first the very simple model function 

AL = asin(V m ) (42) 

without x dependence and the version with \ dependence 

AL = a sin(V>i„)/(l - c{x in ~ Xo)) (43) 

The small perturbative parameters b and c will be considered as of the same 
order. 

Now we consider the preimages of these scattering model functions. First, as 
simplest possibility combine Eqs. [40] and [42] In equation [40| we find 2 possible 
real preimage values of x as l° n g as Pout < Po and real preimage values if 
Pout > Po- Exactly at p out = po two preimages collide and turn from real to 
imaginary, as physical solutions they disappear. Accordingly the line p ou t = Po 



is a caustic line in this case. In equation 42 we find 2 real values of tp for \ AL\ < a 



and zero preimages for |AL| > a. Along the lines AL = a and AL = —a the 
two solutions collide and turn from real to imaginary. Therefore also these 
two lines are caustic lines for this simple case. In total we find 4 preimages 
inside of the rectangle delimited by the lines p out = 0, p out — po, AL — a and 
AL = —a and zero preimages outside of this rectangle (see figure [6]). Of course, 
the caustic structure of this very simple case is not structurally stable, it changes 
qualitatively under small deformations of the scattering functions. Along generic 
caustic lines the number of solutions changes by 2 and not by 4. Therefore we 
need to add the appropriate perturbations to the scattering functions to turn 
the caustic structure into a structurally stable one. We have to see next that the 
transitions from equation [40] to equation [41] and from equation [42] to equation 
|43| are the appropriate perturbations to turn the caustic lines into generic and 
structurally stable ones. 



Next, let us combine Eqs 41 and 42 In this case the function for AL still 



produces the degenerate caustic lines AL = a and AL — — a. If we invert 
equation [41] and we eliminate the ip dependence by inserting from equation [42] 
then we obtain for x the equation 



= xo ± \M> " Pout ± by/1 - AL 2 /a 2 (44) 
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First we see the caustic lines AL — a and AL = — a which are already known 



from the inversion of equation 42 In addition we also see that solutions collide 



and turn from real to complex along the ellipse given by the equation 

a 2 (p Q ~ Pout) 2 + b 2 AL 2 = a 2 b 2 (45) 

In total for the combination of equations [41] and [42] the behaviour of the preim- 
ages is the following: Outside of the strip delimited by the lines AL — a and 
AL = —a the number of real preimages is zero. Let us now look at the interior 
of this strip. For values oip ou t outside the ellipse and at the side of larger values 
we find zero real preimages in x, inside the ellipse we find 2 preimages and for 
values of p out between zero and the ellipse we find 4 preimages. The result is: 
The perturbation introduced in equation [41] is able to modify the degenerate 
caustic line p ou t = Po into a structurally stable curve. However, the two caustic 
curves AL = a and AL = —a stay in their structurally unstable form. 

Next let us combine Eqs(40] and [43j The equatior]40] has the caustic line 
Pout — Po as before. Inserting equatior]40] into equatior]43] leads to the following 
equation for ip 

sin^m) = AL(1 ± cy/p -p ou t)/a, (46) 
which in turn leads to the caustic curves 

AL = ±a{l±c^p a -p out )- 1 (47) 

In the relevant region of p values the curves given by this equation come close 
to parabolas centred at ±a. Now the preimage structure is the following. We 
find zero real preimages for p out > po- Therefore let us concentrate on the strip 
of Pout values between zero and the line p ut = Po- In the region of small values 
of AL between the two approximate parabolas we find 4 real preimages, inside 
of the two approximate parabolas we find 2 real preimages and in the rest of the 
strip zero real preimages. The perturbation introduced in equatior ]43| is able to 
turn the degenerate caustic lines AL — a and AL — —a into structurally stable 
curves. 

Last, let us check that the combination of Eqs 4l] 43 produces structurally 



stable caustic curves only. Eliminating the variable ifji n and making the substi- 
tution x = Xo ~ Xin we g e t the following polynomial in x 

P(x) = x 4 + Ax 2 + Bx + C (48) 

where 

A = 2{ Pout - p ) + b 2 AL 2 c 2 /a 2 (49) 

B = -2b 2 AL 2 c/a 2 (50) 

C = ( Pout -p ) 2 + b 2 (AL 2 /a 2 -l) (51) 

In the following we use the abbreviation D = AL 2 /a 2 — 1. A polynomial has 
colliding solutions whenever its discriminant is zero, see e.g. section II. 4 and in 
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particular lemma 4.7 in [53]. In the case of the polynomial of 48 the discriminant 
is given by 

Dis = -27B 4 - 4A 3 B 2 + 16A 4 C + 144AB 2 C - 128A 2 C 2 + 256C 3 (52) 

We will evaluate the resulting complicated expression in p and x perturbatively 
in the small quantities b and c. The lowest order contributions come in order 4 
and are 

Di Si = 256b 4 { Pout - Po ) 2 D 2 (53) 

This discriminant is zero along the lines p ut = po, AL = a and AL = —a. All 
these 3 lines come with multiplicity 2. It is the same caustic structure as we 
obtained with the combination of the Eqs|4"0"]and[42[ 

Next we keep all contributions up to order 6 combined in b and c and leave 
out the irrelevant constant factor 256& 4 . The result is 

Dis 6 = b 2 D 3 + {p out - Po ) 2 D 2 + 2( Pout - p ) 3 (AL 2 /a 2 + l)AL 2 c 2 /a 2 (54) 

The real zeroes of these last expressions in the (p ou t, AL) plane are sketched in 
figure [6j We have also produced a theoretical plot for the curve of the zeros of 
the expression [54] with the parameter a set on the value 1 whereas the small 
parameters b and c take the value 0.1, see the figure [8(d)] In these plots the 
various regions of the plane are marked by the number of real preimages of the 
original equations[44]and[43] Here all structural instabilities of the caustic curves 
are removed. Note that this figure coincides with the projection singularities 
of the ring shaped mountain shown in figure [7] The projection is highlighted 
on the left side. The parts a, b, c of figure [8] show for comparison the rainbow 
structure in the cross section coming from a single interval of continuity of the 
scattering function for the 3 examples introduced in section 2. Part a belongs 
to the channel with obstacle, part b to the map model and part c to the bottle 
billiard. 

We see: For all the 3 examples the rainbow singularity in the cross section has 
the same qualitative structure as the curves described by [54] This motivates us 



to propose 54 as normal form of the rainbow singularity in the double differential 
cross section for the class of scattering systems considered in this article. 

The shape of the normal form induces a ring shaped mountain. Since a ring 
shaped mountain can be considered half a torus, our projection is half of the 
well known projection singularity pattern of the torus, see fig. 6.13 of [35] . 

For moderate breaking of the symmetry each interval of continuity, which is 
still some strip running around the torus in ip direction, produces in the cross 
section one copy of the typical rainbow singularity. Of course, for each individual 
interval it takes different values of the parameters. The total set of singularities 
seen in the cross section coming from all intervals is expected to consist of a 
superposition of a fractal of shifted and continuously deformed copies of the 
normal form. In principle there should be an infinity of them. Practically we 
can resolve a finite number up to some finite level of hierarchy of the underlying 
fractal. Let us turn next to some numerical examples for the cross section. 
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Figure 6: The diagrams of the regions divided by a rainbow singularity. The 
numbers indicate the quantity of preimages of the region for each point. 




(a) (b) 

Figure 7: On the left side, labelled 7(a)| we represent an annular region which 
spans a quadratic maximum. On the right side, |7(b)| we emphasise the singu- 
lar set of the projection into a two dimensional manifold. The domain of the 
Scattering function is labelled D. The numbers represent the number of inverse 
preimages for the projection. Compare with figures |6| and |8(dy| 
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(c) Rainbow Singularity for the billiard. 



(d) As a comparison, the zeros of the 
polynomial Dis§, equation |54| The pa- 
rameters are a = 1, 6 = c = 0.1 



Figure 8: Different isolated rainbow singularities corresponding in each case to a 
single domain of continuity of the scattering functions, for the different systems 
presented and the solution for the zeros of the Diss polynomial. 



A=0.010, LJ = 1 .0383 



< 




Figure 9: The rainbow singularities for the map system generated by equation 
|T6| The parameters are A — 0.010, L in = 1.0383. We show the structure of the 
projection of the half torus repeated as far as resolution goes. 
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AO.015, LJ = 1 .5575 




To obtain the figures [HJ [9] and [TU] we have done a coarse graining of the 
domain of the cross section which is similar to what any real detector in an 
experiment does. The domain is divided into many small rectangles and the 
counts in any one of these rectangles are registered. Each rectangle can be 
considered one detector channel. High count rates are indicated in the figure 
by a darker colour. Note that there are high count rates also in some detector 
channels which do not contain contributions from the true singularity itself, 
in particular the detector channels near the cusp of the inner singularity line 
show this behaviour. They are the ones roughly connecting the inner cusp with 
the outer singularity line and marked by broken lines in figure |7(b)[ The high 
count rates lie along curves which come close to one ellipse and two parabolas, 
compare with figure [7[ 

After studying the projections of these points we recognise the shape of the 
half toroidal mountains, see figures |9j [10] These figures have been obtained 
by simulating a great quantity of trajectories with initial conditions uniformly 
distributed on the torus. Then we have selected points on the histogram of final 
conditions which show a much higher count than their neighbours, as seen on 
the figures [9] and [TUJ If we compare them with the black outline on the figure 
[7J we can see the projection of the same basic repeated over. For a symmetric 
case we can see the domains of continuity in smoothly changing colours in the 
figure [TT] which partitions the torus into a fractal family of stripes. Each has 
typically an extremal set of measure zero and shows up as a rainbow singularity 
in the cross section. As the symmetry gets broken, the embedding remains 



stable, although the discontinuities bend and may form rings (see figs. 11 and 



12), or even multiply punctured domains. Then we observe distinct shapes 



besides the circular rings which should also represent a complete partition of 
the domain. It should be noted that the smoother parts of the function have 
still a ring shaped domain of continuity, and that they are the main contributors 
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-6.96-6.95-6.94-6.93-S.92-6.91 -6.9 -6.S9-6.SS-6.S7 

Figure 11: The scattering function for the map generated by the function in 
expression [T6| with A = 0.000 and £j„ = 1.0383 on the torus of initial conditions. 
As L is conserved, we only show the p ou t component. 



to the cross section. Therefore, the experimentally detectable signature of the 
partitions of the domain is still a family of ring shaped sets. 

6 Final remarks 

By asymptotic measurements we observe which rainbow contributions are present 
or absent and can conclude which corresponding intervals of continuity in the 
scattering function contribute and finally draw conclusions about the structure 
of the chaotic invariant set in the Poincare map of the system. In particular 
we can change a parameter of the system, for example the symmetry breaking 
parameter A, and follow important changes of the chaotic set, at least in the 
lower levels of the hierarchy of the fractal structure. For further analysis of the 
resulting data two promising possibilities exist. 

First, we can try to construct a symbolic dynamics for the system. For sys- 
tems with two degrees of freedom a description of the development scenario of 
the chaotic set in terms of a development parameter related to an approximate 
symbolic dynamics has been presented for binary and for ternary symmetric 
horseshoes in |13j and [36] respectively. This description was based on a rather 
simple approximation for the symbolic dynamics. In the meantime for chaotic 
sets of two dimensional maps, more sophisticated approximations for the sym- 
bolic dynamics have been developed [37j [38] . It would be worth to generalise all 



27 




-6.96-6.95-6.94-6.93-6.92-6.91 -6.9 -6.S9-6.SS-6.S7 



(a) The p ou t component. 




(b) The AL component. 

Figure 12: The scattering function for the map generated by the function in 
expression[T6|with A — 0.010 and L in — 1.0383 on the torus of initial conditions. 
This gives the rainbow singularities depicted on figure [9] 
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(a) The p ou t component. 




-6.96-6.95-6.94-6.93-6.92-6.91 -6.9 -6.S9-6.SS-6.S7 

(b) The AL component. 

Figure 13: The scattering function for the map generated by the function in 
expression[T6|with A = 0.020 and L in — 1.0383 on the torus of initial conditions. 
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these developments for the four dimensional maps encountered in the present 
article. 

Second, we can measure the weight of the contributions to the structures in 
cross sections and scattering functions coming from the various levels of hierar- 
chy of the underlying fractal. Thereby we extract scaling factors of this fractal. 
The distribution of these scaling factors can be analysed by thermodynamical 
methods to extract the statistical measures of the chaos of the systems. For the 
thermodynamical methods see [39l [40] . For the application of these methods to 
chaotic scattering with two degrees of freedom and its cross section see [41] . 

By knowing which rainbow singularities are present and which ones are miss- 
ing compared to the complete case and by knowing the scaling factors between 
contributions from the various levels of the hierarchy we have knowledge about 
the topology of the horseshoe and about the measures of chaos in two dimen- 
sional cases. 

All our considerations have been restricted so far to the case where only one 
degree of freedom is open and all other degrees of freedom are closed. This sit- 
uation implies the following property of the system: For any value of the total 
energy the closed degrees of freedom can swallow all energy such that for the 
open degree of freedom only energy zero remains. Therefore for any positive 
value of the total energy the motion of the open degree of freedom can come to 
a stop at infinity ( when the potential is attractive without outer barrier ) or at 
the outermost potential barrier ( in cases where it exists, where the potential is 
repulsive for large distances ). In this sense the existence of at least one closed 
degree of freedom implies that for any positive value of the total energy the sys- 
tem sits exactly on a channel threshold. As a consequence we have the dividing 
surfaces of dimension 3 ( in the domain of the Poincare map ) formed by such 
trajectories and the intersection of the stable and unstable dividing surfaces 
forms the chaotic set described before. The elementary intersection set between 
two hyper-surfaces of codimension 1 is a set of codimension 2. If the dividing 
surfaces form fractal folds, then we have a fractal collection of elementary inter- 
section patterns and the complete intersection set has a codimension less than 
2. In our particular case this property guarantees that the chaotic invariant set 
in the Poincare map has a dimension larger than 2. 

The situation is different if no closed degree of freedom exists. Imagine a 
system with only open degrees of freedom and an attractive potential. Then 
for positive values of the total energy it is impossible to have trajectories which 
go out arbitrarily far and where at the same time the velocity goes to zero. 
However, such trajectories are exactly the ones which form the homoclinic tangle 
in the case investigated in the present paper and cause the fractal structure in 
the scattering functions and in the cross section. This consideration shows 
that the structures described in this article will not necessarily exist in systems 
without closed degrees of freedom. 

Of course, topological chaos can exist in systems with only open degrees of 
freedom. One system of this type with 3 degrees of freedom, investigated in 
the past [25l |42] , is scattering of a point particle off four hard spheres situated 
at the four corners of a tetrahedron. When the radius of the spheres is small 
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compared to their distance then all periodic orbits are completely unstable and 
we have hyperbolicity. In the 4 dimensional Poincare map this implies the ex- 
istence of hyperbolic periodic points but excludes the existence of NHIMs. The 
invariant manifolds of the hyperbolic points are 2 dimensional, their elementary 
intersection structures are points. Even when the manifolds are folded to form 
some fractal structure, the complete intersection structure between stable and 
unstable manifolds is a fractal built up of points. Then in general the dimen- 
sion of this intersection structure is small and does not cause much observable 
structure in scattering functions. According to |25j the intersection structure in 
the 4 dimensional domain of the Poincare map should have at least dimension 
2 in order to cause easily observable effects in generic scattering functions. 

The choice of our class of systems has also favourable consequences for the 
cross sections. Naturally we look at the doubly differential inelastic cross section 
as function of two action type variables. In section 5 we studied the cross section 
as function of angular momentum transfer and of outgoing momentum of the 
open degree of freedom. Because of conservation of total energy and because of 
the monotonous dependence of the energy of the oscillating degree of freedom 
on its action we could equally well express the same cross section as function 
of angular momentum transfer and the final action of the oscillating degree of 
freedom. The use of the outgoing momentum has the advantage that its sign 
indicates immediately whether a particular trajectory describes transmission 
or reflection. Therefore we stay with the variable p ou t in the cross section. 
However, because of the above mentioned considerations we treat p ou t as if it 
would be an action variable. Because of energy conservation these two action 
like variables of the cross section are restricted to a finite interval of possible 
values. On the boundaries of intervals of continuity of the scattering function 
Pout goes to zero. In the interior of intervals of continuity the scattering function 
is smooth. Then this function necessarily has generic extrema in the interior 
and has lines along which the determinant of derivatives is zero. This guarantees 
the existence of the rainbow structures described above. 

Again the situation can be different for scattering systems with open degrees 
of freedom only. For simplicity think of the scattering of a points particle off a 
localised potential in a 3 dimensional position space. For the moment consider 
the rotationally symmetric case where the azimuth angle is irrelevant. Then the 
possibility exists that in each interval of continuity of the scattering function the 
deflection angle goes monotonically from minus infinity to plus infinity without 
any generic extremal points. As has been pointed out in |43j this situation also 
can happen in cases of chaotic systems. Then the fractal chaotic set in the phase 
space does not leave fractal traces in the cross section. 

What happens for even more degrees of freedom, let say N? We can make 
some comments for the class of systems where we have one open degree of 
freedom coupled strongly to one closed degree of freedom and any number of 
further degrees of freedom which are only coupled weakly to the first two de- 
grees of freedom. Let us assume again that there is a parameter A which gives 
the strength of this coupling and that we have again a limiting case A = 
where these additional degrees of freedom are decoupled from the first two ones 
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and for simplicity also decoupled among themselves and where accordingly the 
system can be reduced to one with two degrees of freedom. In this case we find 
a (2N-4) dimensional NHIM in the (2N-2) dimensional domain of the Poincarc 
map. It is stable in one direction, unstable in one direction and neutrally stable 
in the remaining (2N-4) directions. The development degree of the horseshoe 
of the reduced system depends on the amount of energy in the reduced sys- 
tem. Since it depends on the value of the conserved actions of the remaining 
degrees of freedom we can use one of the remaining neutral directions as de- 
velopment parameter of the horseshoe and in total obtain again a pile of two 
dimensional horseshoes similar to the 3 dimensional tangle described before for 
the three degrees of freedom case. In particular we can apply the same argu- 
ments of robustness as before. The difference is that now we have to form the 
Cartesian product of this three dimensional structure with (2N-5) neutral direc- 
tions, which are partially action directions and partially angle directions. The 
domain of the scattering function is now a (N-l) dimensional torus of initial 
relative phase shifts between the N degrees of freedom and its range is a (N-l) 
dimensional interval of final actions. Take our previous quantities L, 9, \ an d 
AL to be (N-2) component quantities. Under small breaking of the symmetric 
case A = the robustness argument indicates again that the coarse structure 
of the scattering function should be stable, the lower the level of hierarchy in 
the fractal structure, the more stable the situation should be. The (N-l)-fold 
differential cross section is the relative probability to obtain some combination 
of final actions for a constant density of the initial phase shifts. The rainbow 
singularities are again the projection singularities of the graph of the scatter- 
ing function under projection on the range. The elementary rainbow structure 
coming from a single interval of continuity of the scattering function is now a 
(N-2) dimensional surface in the (N-l) dimensional domain of the cross section. 
We plan a more detailed discussion of the general case in a future publication. 

So far everything has been explained for classical dynamics. Therefore the 
question remains how our results are reflected in quantum systems, like the ones 
presented in [33]. We have analysed the cross section as function of the out- 
going action for fixed total energy. This was appropriate since classically the 
actions of the closed degrees of freedom are continuous variables. In this point 
quantum mechanics is essentially different. In our models asymptotically the 
open degree of freedom is decoupled from the two closed degrees of freedom, see 
for example equation [2] Accordingly in the asymptotic region the transverse 
motion, i.e. the state of the two closed degrees of freedom, must be in one of 
the discrete quantum states of this bound subsystem. For example in equation 
[2] the transverse motion is a two dimensional harmonic oscillator with its usual 
quantization of the action according to / = (n + l/2)hu. For other particular 
models of the closed transverse degrees of freedom similar restrictions apply. 
On the other hand this discreteness of the asymptotic transverse states provide 
a natural channel decomposition of the S-matrix and the scattering amplitude 
where we start from a particular asymptotic initial transverse state and study 
the transitions to the various other energetically allowed final asymptotic trans- 
verse states. The open longitudinal degree of freedom does not have similar 
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restrictions, its asymptotic energy can be any positive value. Thus to scan the 
action of the closed degrees of freedom is impossible in quantum dynamics. The 
appropriate procedure is to select particular initial and final states of the closed 
degrees of freedom and to vary the total energy continuously, i.e. study the 
cross sections for the various channel to channel transitions as function of the 
total energy. In these quantities we expect to see structures, in particular se- 
quences of scattering resonances, related to the classical chaotic set. So far we 
did not yet study such implications and they are not the subject of the present 
publication. 

The essential point for the analysis of the quantum mechanical cross section 
is the wave dynamical analog to classical rainbows. One essential difference be- 
tween classical dynamics and quantum dynamics is how the contributions from 
the various initial conditions are summed. In classical dynamics the cross sec- 
tion of equation [39] itself is a sum over the contributing initial conditions. In 
a semi-classical approximation to quantum dynamics or when using Feynman's 
path integral version of the quantum propagator, we first form a scattering am- 
plitude as a sum over contributing trajectories and then form the cross section 
or scattering probability as absolute square of the amplitude. Therefore the 
resulting cross section is a double sum over the trajectories and contains first 
a sum over diagonal terms which has the structure of the classical cross sec- 
tion plus the non diagonal terms which are interferences between the various 
contributing trajectories. If in a rainbow two trajectories fall together then in 
a semi classical treatment of the scattering amplitude we have to apply some 
uniformization which replaces the contribution of the two trajectories by an 
Airy function. This gives the wave dynamical analog of the classical square root 
singularity. 

Accordingly two things might be done in a wave dynamical treatment. First, 
the interference pattern in the wave cross section can be analysed. For some 
attempts in this direction see [41] and [45], where the semi-classical scattering 
from three soft potential mountains has been used as example of demonstration. 
Second, a decomposition of the scattering amplitude into Airy type contribu- 
tions might be tried, in order to recover at least the first few hierarchical levels 
of the fractal pattern of rainbows. To our knowledge this has not yet been done 
so far, but might be worth to try. 
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